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Abstract 

We find all polynomials Z{z) such that the differential equation 



in 

|x(z)^ + Y{z)f^ + Z{z)^S{z) = 0, 

Qh. where X{z),Y{z), Z{z) are polynomials of degree at most 4, 3, 2 respectively, has 

^ . polynomial solutions S{z) = HILil^ — ^i) of degree n with distinct roots Zj. We 

-4— > . 
C^ . derive a set of n algebraic equations which determine these roots. We also find all 

polynomials Z{z) which give polynomial solutions to the differential equation when 

the coefficients of X{z) and Y{z) are algebraically dependent. As applications to 

h>. ■ our general results, we obtain the exact (closed-form) solutions of the Schrodinger 

^^ ■ type differential equations describing: 1) Two Coulomb ically repelling electrons 

^^ ■ on a sphere; 2) Schrodinger equation from kink stability analysis of (j)^-type field 

• ' theory; 3) Static perturbations for the non-extremal Reissner-Nordstrom solution; 

(^ • 4) Planar Dirac electron in Coulomb and magnetic fields; and 5) 0{N) invariant 

decatic anharmonic oscillator. 
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1 Introduction and main results 

Consider the general 2nd order linear ordinary differential equation (ODE) 

X{z)^^+Y{z)^^+ Z{z)^ S{z) = 0, (1.1) 

where X{z), Y{z), Z{z) are polynomials of degree at most 4, 3, 2 respectively. 



X{z) = J2 <^kz\ Y{z) = Y, hz\ Z{z) = J2 

k=0 k=0 k=0 

1 



CfcZ^ 



The ODE fll.ip has as many as 12 parameters a^, 6^, c^, and contains, as particular cases, 
the Heun and generahzed Heun equations as well as various confluent equations. As 
examples, here we list five ODEs for each of which we provide a physical application in 
section |3l 

(i) 04 = 63 = C2 = 0, i.e. degX(2;) = 3, degY{z) < 2 and degZ{z) < 1. If X{z) has no 
multiple roots, we can write X{z) = Y[s=ii^ ~ '^s), j^ = X]s=i l^d" ^^r suitable 
complex numbers ds,as- We then obtain the Heun equation. 



(ii) degX(2;) = 4, degF(2;) < 3 and degZ{z) < 2. li X{z) has no multiple roots, we can 



write X{z) = Yl^^iiz — e^), ^^ = ^^^^ ^3^ for suitable c-numbers es,/^^- Then 



(11. ip takes the form. 



^ + V^^:^+ /^^^ \siz) = o. (1.3) 



This is the generalized Heun equation, 
(iii) degX(2;) = 3 (i.e. 04 = 0), degY{z) < 3 and degZ{z) < 2. If X{z) has no multiple 

Y{z) ^ Y"3 _us_ 
X{z) l^s=l z-fs 



roots, we can write X{z) = ns=i(^ "~ /s); yu) ~ Ss=i 'l~r ~'~ ^ ^°^ suitable 



c-numbers fsiVsi ^- Then (II. ip has the form, 

S{z) = 0. (1.4) 




' ^ ' Ziz) 




dz' \tiz-fs dz nLi(^-/. 



This is the equation considered by Schafke and Schmidt [T]. 
(iv) degX(;z) = 2 (i.e. 04 = 03 = 0), degy(z) < 3 and degZ{z) < 2. If X{z) has no 
multiple roots, we can write X{z) = (2; — gi){z — (72), yu) ~ '^ — '' '^ — \- az + k 
for suitable c-numbers gi,g2,cri,a2, a, n. Then (II. ip is given by 

— + —^ + —^ + az + n]- + \^ ^ S{z) = 0. 1.5 

dz^ \z-gi z-g2 J dz [z - gi){z - g2) } 

(v) degX(2) = 1 (i.e. a^ = a^ = a2 = 0), degY{z) < 3 and degZ{z) < 2. Write 
— ^ — ^^ + Xz'^ + 7z + (5 for suitable c-numbers h, rj, A, 7, 6. Then (II. ip reads 



X(z) 



;^+(7^ + ^--^ + -'-- + ^)s+(|^^5W-°- (!■«) 



Recently there is a lot of research interest in finding polynomial solutions to second 
order differential equations of the form (II. ip [2]- [H]. ODEs with polynomial solutions are 
often called quasi-exactly solvable and have wide-spread applications in physics, chemistry 



and engineering (see e.g. [12] and references therein). One of the classical problems about 
the ODE ([n]), suggested by E. Heine [13], T. Stieltjes [H] and G. Szego [I5], isQ 

Problem Given a pair of polynomials X{z), Y{z) and a positive integer n, 

(a) find all polynomials Z{z) such that the ODE U.l\) has a polynomial solution S{z) 

of degree n . 
(h) findS{z). 

In this paper we solve this problem by means of the so-called Functional (or Analytic) 
Bethe Ansatz method [T614T9] . Precisely, we find the explicit values of the coefficients 
C2, Ci, Co of Z{z) which give rise to degree n polynomial solutions S{z) of 01. ip with roots 
Zi, 2:2, ■ ■ ■ , -Zn of multiplicity one, and we obtain a set of n algebraic equations (the so-called 
Bethe ansatz equations) which determine these roots.. 

For cases (i) and (ii) above with fixed real ds-, e^ and positive a^, /i^ numbers in (11.21) and 
(II. 3p . respectively, there is a classical result, known as Heine- Stieltjes theorem [13]. This 
theorem says [5] that if the coefficients of X{z) and Y{z) are algebraically independent, 
i.e. do not satisfy any algebraic relations with integer coefficients, then for an arbitrary 

(ji _|_ dQgXiz) — 2 \ 
polynomials Z{z) of degree 
/ 
exactly (degX(z) — 2) such that the ODE has a degree n polynomial solution S{z). 

However, even for these two cases, no results about the values of the coefficients C2, Ci, Cq 

of Z{z) seem to be previously known. Furthermore, one may ask 

Question // the coefficients of X{z) and Y{z) do satisfy some algebraic relations with 
integer coefficients, i.e. are algebraically dependent, then how many polynomials Z{z) are 
there which lead to degree n polynomial solutions of the ODE U.l\) ? 

To my knowledge, this question was not answered by Heine and Stieltjes in their 
theorem. In this paper we will also provide an answer to this question as a by-product of 



our general result (theorem II. II below). 

We now state one of our main results of this paper. 

Theorem 1.1 Given a pair of polynomials X{z) and Y{z), then the values of the coef- 
ficients C2,Ci,Co of polynomial Z{z) such that the differential equation U.l\) has degree n 
polynomial solution 

n 

S{z) = l[{z-z,) (1.7) 

with distinct roots zi, ^2, ■ ' ' ? ^n ^'^c given by 

C2 = —n{n — 1)04 — nbs, (1.8) 



* Authors in these references only considered the cases corresponding to (jl.2|) and ()1.3|) . while the 
general ODE (jl.ip here also contains many other cases, e.g. differential equations (|1.4p - (|1.6|) . 



Cl 



- [2(77. - 1)04 + 63] '^Zi- n{n - 1)03 - nb2, (1.9) 



n 



Co = - [2(n - 1)04 + h] ^ -2^^ - 2a4 ^ -Zj-Zj 



— [2(r2 — 1)03 + 62] 2. "^^ ^ """(^ ~ ^)(^2 — iT-bi, (1-10) 

where the roots zi,Z2,- ■ ■ ,Zn satisfy the Bethe ansatz equations, 

> + 4— 3— 2— ^ = 0' « = l,2,---,ri. (1.11) 

]^i Jill 

The above equations / (J.<§] )-/ I7T77]) g'we a// polynomials Z{z) such that the ODE U.l\) has 
polynomial solution ( [i.?] ). 

Remark 1.2 T. Stieltjes gave in [14^ a set of equations satisfied by the roots of S{z) for 
fixed real {dg}, {cg} and positive {as}, {fis} numbers in ( fi.^j) and /il.3\) . respectively. Our 
Bethe ansatz equations U.ll\) reduce to those of Stieltjes under these conditions. However, 
even for these particular cases, the expressions for the coefficients of Z{z), U.8\) - I[TT^) . 
had not been obtained previously. Our theorem li.il determines the explicit form of all 
polynomials Z{z) which give polynomial solution to the general ODE U.l\) . From the 
explicit expressions of C2,Ci,Co above, the number of Z{z) is given by the number of the 
solution sets of the Bethe ansatz equations U.ll\) . 

From (ILSp -f fTTTTj) . we also have, as special cases of the theorem II. H 

Corollary 1.3 (a) If for an arbitrary integer n the coefficients a^ and 63 in the ODE 
( li.ij) are algebraically dependent, 

2(^-1)04 + 63 = 0, (1.12) 

then there are n + 1 polynomials Z{z) with coefficients C2 = n{n — 1)04, Ci = 
-n[{n - 1)03 + 62] IMl^ and 

n n 

Co = -n [{n - l)a2 + 26i)] - 2a4 ^ ZiZj - [2{n - 1)03 + 62] ^ Zi, 

i<j i=l 

such that U.l\) has degree n polynomial solution S{z) ( [i.?] ), where the roots z^ satisfy 
KTWl with 63 = -2{n - 1)04. 
(b) If a4^ = b-^ = C2 = and if for an arbitrary integer n the coefficients a^ and 62 are 
algebraically dependent, 

2(n- 1)03 + 62 = 0, (1.13) 

then there is 1 polynomial Z{z) = —n[{n — 1)03 + b2]z — n{n — l)a2 — nhi such that 
the corresponding ODE has degree n polynomial solution S{z) ( [i.Tp with the roots 
Zi determined by lil.ll]) with 04 = 63 = 0, 62 = —2{n — 1)03. 



Other main results of this paper are the exact (closed-form) solutions of five phys- 
ical systems (see section |3] below), providing interesting examples and applications of 
the differential equations (ll.2p - (ll.6l) . respectively. The existence of such exact solutions 
demonstrates that these systems are quasi-exactly solvable. 

In the appendix we present some explicit formulas corresponding to the differential 
equations flT:2|) - flL6|l . 

2 The Proofs 

Proof of Theorem \1.1\ We prove theorem 11.11 by using the Functional Bethe ansatz 
method. This method was used by us to obtain exact polynomial solutions of certain 
higher order ODEs arising from the nonlinear optical and spin-boson systems [T5|[TU] . Let 

n 

S{z) = \[{z - z.-^ , (2.1) 

be a degree n polynomial with undetermined, distinct roots zi, Z2., ■ ■ ■ ■, z^. We will be 
looking for the values of the coefficients C2, Ci, Cq of Z{z) and the roots zi^ 1 < ^ < "n., such 
that (12. ip is a solution of (II. ID with fixed coefficients of X{z) and Y{z). Substituting 
Z{z) into the ODE and dividing on both sides by Z{z) gives rise to 



n n 

- Co = ia^z'^ + a^z^ + a2Z^ + aiz + oq) } > 

\ ' ^-^ Z — Zi ^-^ 7- — 7- 

1=1 J 

n ^ 

[h^z^ + h^z^ + hxz + 6o) ^ -^^ + C2Z^ + 



■ -, Z Zi Zi Zj 

1=1 J^l ■' 



CiZ 
. , Z- Zi 

The left hand side is a constant and the right hand side is meromorphic function with 
simple poles at z = Zi and singularity at z = oo. The residues of — cq at the simple poles 

z = Zi are 

n 

Res{-co)z=zi = {a4zf + aszf + a2zf + aiZi + ao) ^ — - — 

+b3zf + b2Z^ + biz, + bo. (2.2) 



It can then be shown that 
Res(-co)2=zi 



co-Y^ ^_°^"-^' = J2[a4iz^ + zf + z^Zi + zz^) + a3{z'^ + z^ + zZi) 

i=i ^ ^' i=i 



z — z. 

1=1 

n 



■^a2{z -F Zi) -F ai] ^ 

. , . Zi Zj 



+ Y^ [^3(2;^ + zj + zzi) + b2{z + Zi) + 61] + C2Z^ + ci; 



j=i 



[n{ 


n — 1)04 + nh'i + C2] z^ 


+ 


(2(n - 1)04 + 63) ^ Zi + n{n - l)a^ + n62 + Ci 


n n 

+ {2{n - l)a4 + 63) Y^ z] + 2a4 ^ ZiZj 


n 



+ (2(n - 1)03 + 62) ^ Zi + n(n - 1)02 + nhi, 



i=l 



(2.3) 



where we have used identities, 

n n ^ n n 



n n 9 

Z^ 






-nfn — 11 



n n o 



EE-4t = («-i)5:^<. EE-e 



m 



i=i jyi 



Z^ Zj 



i=l 



i=l jyj 



Zi Z^ 



i)En' + E 

i=l j<j 



Zj Zj . 



That is 



- Co = \n{n - l)a4 + 7263 + C2] 2:^ 



+ 



(2(n - 1)04 + 63) ^ 2;j + ?7.(n - 1)03 + nh^ + Ci 



i=\ 



^ Res (-Co 



i=l 



2;- Zi 



Z^Zj 



+ (2(r2 - l)a4 + 63) XI ^*^ + 2a4 X , 

n 

+ {2{n - 1)03 + 62) X -Zj + ^('"' - 1)«2 + ''^^i- 



(2.4) 



i=l 



The right hand side of (12.41) is a constant if and only if the coefficients of z^ and z as well 
as all the residues at the simple poles are equal to zero, respectively. This gives C2, ci in 
terms of the fixed coefficients of X(z), Y{z) and the roots of S{z), 

n{n — 1)04 + nbs + C2 = 0, 

n 

(2(n - 1)04 + 63) X Zi + n{n - 1)03 + 77,62 + Ci = 0, 

i=l 

and the n algebraic equations determining the roots Zj, 



Y^ 2 63-2/ + ^22;^ + &i^i + bo 



j¥=i 



Zi - Zj a^zf + 032;^* + a2zf + aiZi + ao 



^ = 1,2, 



, ra. 



(2.5) 



These equations are nothing but (11. 8p . fll.9p and fll.lip . respectively. It then follows from 
(El]) that 

n n 

- Co = (2(n - l)a4 + 63) ^ z^ + 2a4 ^ ZiZj 

n 



+ (2{n - 1)03 + 62) 2J -^* + ^(^ ~ 1)'^2 + "-&1 



=1 



which is (ll.lOp . The proof is thus completed. 



D 



Proof of Corollary \1.3[ Part (b) is obvious. We now prove part (a) . If the coefficients of 
X{z) and Y{z) are algebraically dependent, 63 = —2{n — 1)04, then the ODE (11. ip has a 
hidden sl{2) algebra symmetry if C2 = n{n — 1)04, C\ = —n[{n — 1)03 + 62]- To show this, 
we write the ODE (II. ip with such coefficients C2, Ci as the Schrodinger form 

HS{z) = -coS{z). (2.6) 

Then it can be shown that H is an element of the enveloping algebra of Lie-algebra s/(2), 

H = a^J^ J^ + a^J^ J"^ + a2J^ J"^ + aiJ^ J- + a^J- r 
'I 



+ 



^,(3^_2)a3 + &2 



where 



J+ + [{n - 1)02 + 61] J° 

(Ti \ ft ft 

-ai + 6oj J' - -^ci2 + 2 [(^ - 1)«2 + hi] , (2.7) 

J^ = z'^-nz, J' = z^-^, J- = ^ (2.8) 

az dz 1 az 

are differential operator realization of the n + 1 dimensional representation of the s/(2) 
algebra. — Cq is the eigenvalue of H with polynomial eigenfunction S{z) = YYi=ii^ ~ ^i)i 
given by (from theorem II. ip 

n n 

- Co = n [(n - l)a2 + 2hi)] + 2a4 J^ -Zi^;^ + [2(n - 1)03 + 62] ^ ^i, (2.9) 

and the roots Zi of the eigenfunction are determined by the Bethe ansatz equations, 

i:^+ """:''°fi^''ri^''i^''° =o. .=1.2....... (2,10) 

It is well known that the solution space of second order differential operators with a s/(2) 
algebraization is n + 1 dimensional [8l[9]. Applying this to the H above, we conclude 
that the above Bethe ansatz equations have n + 1 sets of solutions and thus there n + 1 
eigenvalues — co, i.e. n+1 polynomials Z{z). 



D 

Let us remark that as far as we know the Schrodinger equation fl2.6p with H given by 
(12. 7p had not been solved previously (except for some very special cases), and thus (12.91) 
and fl2.10p above give the first exact solution for the general ODE (12. 6p . 

3 Applications 

In this section we apply our general results obtained in section 1 to derive exact solutions 
of five physical systems, providing examples and applications of the differential equations 
(ll.2p - ll.6l) . respectively. 

3.1 Two Coulombically repelling electrons on a sphere 

Consider a system of two electrons, interacting via a Coulomb potential, but constrained 
to remain on the surface of a D-dimensional sphere of radius R [20]. The Hamiltonian of 
the system (in atomic units) is 

^ = 4(^1 + ^')-^' (3.1) 

where u = |ri — r2| is the inter-electronic distance. The Schrodinger wave function of the 
system can be separated as a product of spin, angular and inter-electron wave functions. 
The inter-electron wave function '^(u) satisfies the ODE 



-hr- + - = ^^' (3-2) 



where 8 and 7 are parameters related to the dimension D of the sphere. Introduce 
dimensionless variable z = -^. Then the above ODE can be written as 

d?^[~^ z+1 + z-1 J d~z + z{z + l)(z-l)f ^ - "• ("^-"^^ 

This ODE has the form of the Heun equation (11.20 . It follows from our general results in 
previous sections that this equation has polynomial solutions of degree n = 1, 2, • • • , 



n 



^{z) = l[{z-z,), (3.4) 



i=l 



where Zi are the roots of the above polynomial to be determined, provided that E and R 
take the values given by 

Tl 

E = -^{n + S-1), (3.5) 

1 " 

R = --[2{n-l) + 6]Y,^^, (3.6) 

8 



where zi, Z2,- ■ ■ ,Zn obey the Bethe ansatz equations 



y2 2 , 1/7 , U^-ih) , U^-ih) 



Zi + l 



Zi - 1 



0, i = 1, 2, ■ ■ ■ , ra. 



(3.7) 



The energy eigenvalues E agree with those obtained in [20] by a different method. How- 
ever, the above general, exact formula for the radius R had not been given previously. 
For n = 1, we have, from the Bethe ansatz equations, Zi = ±^=- Then 2R = —6zi = 

\/~ (by choosing the negative root zi = — -^ so that the radius R is non-negative), 
E = 'J and the wave function is \l/ = z + -h= = ^=(1 + 'ju). 
For n = 2, we find 



Zl 



Z2 



2(, + 2),-V2(^ + 2) + ^^V'(^ + 2)-' 



1 



2(5 + 2) 



' 2{b + 2) + ^^^^ T j2i6 + 2)-- 



(3.1 



The other root satisfies zi = —Z2 which leads to -R = and therefore is discarded. The 
radius R and the energy E are 



R 

E 



-\{6 + 2){z, + Z2) = \j2i6 + 2) + ^^, 
7(5 + 1) 



7(5 + 2) + 25 + 3 



and the wave function is 

^ = {Z-Zi){z-Z2) 



7(5 + 2) 



1 + 7U + 



7^(5 + 2) 



27(5 + 2) + 45 + 6 



u 



(3.9) 



(3.10) 



The ODE (13. 3p has a hidden s/(2) algebra symmetry. To see this, we rewrite (13. 3p as 

H^ = -2R^, H = z{z^ - 1)4^ + (5^2 - 1) A - AR^Ez. (3.11) 

dz'^ 7 dz 

If E takes values given by (13. 5p . i.e. AR^E = n{n — 1 + 5), then H is combination of sl{2) 
algebra generators ( 12. Sp . 



H = rr - rr + 



-(3n-2) + 5 



J^ 



1 '^'i T- 



(3.12) 



This provides an sl{2) algebraization of the two electron system. This algebraization had 
not been realized previously. The eigenvalues —2R of H with the polynomial eigenfunction 
\& = niLil-^ ~ ^i) ^^^ given (13. 6 p and the roots Zi are determined by the Bethe ansatz 
equations (13. 7p . 



9 



3.2 Schrodinger equation from the kink stability analysis of the 
0^-type field theory 

Consider the 0^-type field theory in 1+1 dimensions characterized by the Lagrangian 



C^\d., 



f^^ 



2j,2\2 



B^2(l + e2 



(/0^ + e') (1 - /0^) 



(3.13) 



where e is real dimensionless constant and // has the dimension of mass. As is well-known, 
the field equation of this theory has a kink solution. When one performs stability analysis 
around the kink solution one arrives at the Schrodinger equation 



d^ 



dx"^ 



+ V{x) \ip{x) = Eip{x), 



(3.14) 



where E > and the potential V{x) is given by 

,8sinh^ f - (f - 4) sinh^ £^ + 2 (^ + l) (^ - 2) 



V{x) = /i' 



8(l + i + sinh2f)' 



(3.15) 



We show that the Schrodinger equation can be transformed into the form of (II. 3p . To 
this end, we let 






1 + — + sinh^ — 
2 



y 



Then it can be shown that y satisfies the ODE 
sinh^cosh^ ^^-^ 



y +'^^^— ^^^+IT4 + sinh 



l + ^ + sinh'^' 



f^' (1 + ^) 

2 fixy 



E 



r/^' ) y- 



Make a change of variable. 



Then the above equation becomes 



-*+l^-i)-'^-^ 



z = cosh 



fix 



(3.16) 



(3.17) 



(3.18) 



1 
AE 



5z-'- {- + 6]z 



I 






^-6 



y = 0. 



(3.19) 



This ODE can be easily cast into the form of the generalized Heun equation (II. 3p . noting 
that X{z) = z^ + {^ - l) z^ - J2 = {z + l){z - l){z + ^)(z - ^) has no multiple roots. 
From our general results in previous sections, the ODE (13.19^ has polynomial solutions 
of degree n = 1, 2 ■ ■ ■ , 

n 

y{z) = l[{z-z,), (3.20) 



j=i 



10 



where Zi are the roots of the above polynomial to be determined, provided that £", e satisfy 
the relations 

,2 



E 



/^ 



{n — 1)(5 — ra), 



i=l 

6(n-l) 



\n 



l)(n - 6) + [5 -2(n- 1)1^^,2-2^ 



Z-i Zj , 



(3.21) 
(3.22) 

(3.23) 



1=1 



i<j 



and the Bethe ansatz equations, 



4 + 



I z - 



E 



5zf 



j^i * ^ 



:^ + 6] Zi, I 



1,2, 



n. 



For n = 1, we have ii^ = and zi = so that y{z) 



and 



ip{x) 



coshf 



(l + i + sinh^f)^ 
which is the ground state of the system with energy eigenvalue E = {] 

3 

4 



(3.24) 

There is no constraint on e, 

(3.25) 



For ra = 2, we have E = |/x^ and equations. 



-^ = -4 + 3{zl + zj) - 2z^Z2, 



^i + >,2 



4+>^. 



lU - 



1 62 



1 Uo^ - 



' 62 



2;i -^2 
2 

2:2-^1 



Zl 



5zl 



hzi 



-Z2, 



+ 6U1 



+ 6^22 



(3.26) 
(3.27) 

(3.28) 



which give zi = —Z2 = \/2 and ^ = 2. Thus y{z) = (z — \/2){z + -\/2) = z'^ — 2 and 

cosh^f -2 ^ ^ 

ijix) = ^ ^. (3.29) 



(l + ^ + sinh^f) 



3 ■ 

2 



This gives the first excited state eigenfunction. The two energy eigenvalues and eigen- 
functions above for n = 1,2 reproduce those obtained in [21] (from a completely different 
analysis). Here, we have obtained the closed form expressions for all eigenvalues and 
eigenfunctions of the system. From (I3.2ip . we see that there are only five non- negative 
energy solutions to the Schrodinger equation with energy eigenvalues E < fi^. All other 
analytic solutions (corresponding to ra > 6) have negative energy eigenvalues. These 
negative energy solutions correspond to unstable modes. 



11 



3.3 Static perturbations for the non-extremal Reissner-Nordstrom 
solution 

Consider the Lagrangian describing 4- dimensional gravity coupled to an abelian gauge 
field Afj^ and a real massive scalar (with mass rris) [23] , 



C = V9 



mG 2^"^^^^ 4 

1 o ,o „,,. 1 






WtcG 2' '^^' 4 

V{<jy) = ^m^02, /(</.) = -^— ., (3.30) 



where R is the scalar curvature, G is the Newton constant, and a is a parameter which 
has the dimension of length. The interaction between the gauge and scalar fields is non- 
renormalizable. For static solution with magnetic charge Qm, defined by J^2 F = 4:7Tgm, 
where F = ^F^^dx^ A dx'^ is the gauge field 2-form, one must have [23] 

ds^ = gttdt^ + Qrrdr'^ + r'^{d9^ + sin^ 9 d(p^), 

F = g^de A sine d<^, = 0(r). (3.31) 

The relevant equations of motion are 

G^, = 8nGT^,. (3.32) 



By expressing g^ = —e^^^'^'^ and grr = e'^^^'^\ one obtains (in the unit I/vSttG = 1) the 
following system of equations for (p and B, 



0'2 ^ ^2By^^ _ _^/ ^ ^ ^ Q ^333) 



1 2 1 - p2B 

U'' + e^^Kff -iB' + ^ 

Linearizing the above equations around = and fitting e^^^ to the Reissner-Nordstrom 
form for large r gives rise to [23j 

^"+0_2B')^' = e-^(O,r)^, 

ocl + ^^^^^^;^ -(r - 1 for r ^ 1, 

^m-2 

0oce-™^" for r->oo, (3.34) 

where M is the mass of the Reissner-Nordstrom black hole. In the non-extremal case 
M^ > 327r^5f^, the function r'^ - ^r + ^ has two roots 

r± = ^±^VM2-327r2^^, (3.35) 
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p{r) = 


r-r+ 


+ - 
r 


7 

— r_ 


q{r) = 


2 1 

= -m^ + 


2a2 
r2 


+ 2a2 



which correspond to the two horizons of the black hole. Noting that r± obey the relations, 

r+r- = ^, r++r^ = ^, (3.36) 

then the ODE for can be brought into the form, 

0" + p(r)0' + g(r)0 = O. (3.37) 

Here 



•J_ + ±)l + ^^ + ^?_ (3.38) 

'' ly ly iyt ly ty ry ry 

with 

q- = -^^(my^ + ^^ + ^{l-r+ + r_)]. (3.39) 

The Lagrangian has a rigid rescaling symmetry. This scaling freedom can be used to set 
the horizon location r+ = 1. Then black hole mass M can be determined by requiring 
that the horizon is at r+ = 1. It follows from (I3.36P that M = 27r {g"^ + 2) and r_ = ^. 
The corresponding g(r) and p{r) reduce to HI 

P{^) = 7 + 



1 ' 

qrt I jj-i jj-i 

2a^ 1 1 

g(r) = -m^ + ^- + 2a^(l + — )- 

T^ Z /y /y 



2 



a'^^-^Kl + r!) 1 ^ ^'s9lr^ + V7 1 ^3^^^^ 



n fy ry I T ry ry /y 

We now find exact solutions to (I3.37P with g(r), p{r) given by (13.401) . By means of 
transformation 

0(r) = r" e""'^" /(r), ^ = ^ (l ± Vl - 8aA , (3.41) 

then it can be shown that /(r) satisfies the following ODE 

/ + — + 7 + 2ms] f + -, -r, r/ = 0, 3.42 



t Authors in [25 also considered the differential equation (I3.37P but arrived at the different values of 
M, r_ as well as the different p(?'), <?(?') functions. 
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where 



1 1 / 2a^ 

C2 = 2a2(l + — )-2m,(/i + l) + - [ g'^^{o? + mlr _) - ml{l + rV) + 

r_ 1 — r_ \ r_ 

ci = K(2/i + l) + ^](r_ + l) "— '- 

g^{a +mjr_ -m^(l + r_)r_ + 



Co = 2a2 + 2 [a2-^(m, + 1)] r_. (3.43) 

This ODE is of the form (11 ■4p . Applying our general results in previous sections, we can 
show that ( 13.42P has polynomial solutions of degree n = 0,1,2,- ■■ , 



n 



fir) = JJ(r - n), / = 1 forn = 0, (3.44) 

where r^ are the roots of the above polynomial to be determined, provided that a, m^, Qm 
obey the following relations: 

2a2(l + — )-2m,(/i+l) 
r_ 

1 / 2a^\ 

gl^{a^ + mlr^) - m^(l + r^) + = 2m,n, (3.45) 



l-r_ V^"^ ■ . / «v ■ -/ ■ ^ 

2a2(r_ + 1)2 



K(2/i+l)+^](r_ + l) 



r_ 



1 / 2a2 

Qmia^ + m1)r_ - m^(l + rty^ + 



n 

= 2m, ^ r, - n[n + 2/i + 1 + 2m,(r_ + 1)], (3.46) 

2a^ + 2[a^ - fi{ms + 1)] r_ 

n n 

= 2ms ^ r^^ - 2 [n + /i + ms(r_ + 1)] ^ 



n 

+ n[{n- 2/i) (r_ + 1) - 2m,r_] , (3.47) 



where the roots Tj are determined by the set of Bethe ansatz equations, 

V^— + ^ + ^- + — i 2m, = 0, i = l,2,---,n. (3.48) 

^ — ^ Tj — Tj Ti Ti — 1 Ti — r_ 

Note that /(r) = 1 is a solution of the ODE (I3.42p provided that a, rus, Qm satisfy the 
following relations 

2a\l + —)- 2m,(/i + 1) + -^- ( gl{a' + mlr^) ~ ml{l + r! ) + ^) = 0, 
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K(2/i + l)+li] (r_ + 1) ^—^-L 



gl{a^ + my.-ml{l + rl)r^ + — ]=0, 



l-r_ V .V . -. . ^_ 

2a^ + 2[a^- fx{m, + 1)] r_ = 0. (3.49) 

These relations can be obtained from (I3.45p - (l3.47p by letting n = 0. It follows that 

0(r) = r'2(^±VT^8^)^-msr^ (3^5Q) 

where a and m^ are determined by equations fl3.49p . This gives the first exact solution of 
the differential equation (13.371) such that the horizon is at r_|_ = 1. 

3.4 Planar Dirac electron in Coulomb and magnetic fields 

Consider the (2+l)-dimensional relativistic system of a Dirac electron (with mass rrie) 
in the presence of an external electromagnetic field A^. This system was also examined 
in [25] via a similar Bethe ansatz approach. The covariant Dirac equation (in the unit 
h = c = 1) has the form 

t^id^ + teA^)^{t, r) = m.^it, r), (3.51) 

where rrie is the rest mass of the electron, — e (e > 0) is its electric charge and the 
(2 + 1) Dirac gamma matrices 7^^ satisfy the anti-commutation relations {7^,7*^} = 27]^" 
with 1]'^'^ = rj^i, = diag(l, —1, —1). In an external Coulomb and a constant homogeneous 
magnetic field B, the vector potential can be written as 

Ze . By , Bx 

Ao = , ^1 = -^, ^2 = —. 3.52 

r 2 2 

Then the Hamiltonian H{r) of the system can be expressed as 

ido-^{t, r) = i7(r)^(t, r), H{r) = 7V^fc + eAo + 7 V, (3.53) 

where Pk = —idk + eAj., k = 1,2, is the operator of generalized momentum of the electron. 
The wave function \l/(t, r) is assumed to have the form 

^{t, r) = ^ exp{-iEt) iJi{r, 9), (3.54) 



where E is the energy eigenvalue of Hamiltonian, and 

F(r) e*'^ 
G(r) e^('+i)^ 



F(r) e*'^ 
Mr,e)=\ _^„,,, 1, (3.55) 
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where / is an integer. Substituting fl3.54p and f l3.55p into f l3.53p and working in the polar 
coordinates (t, r, 6) reduce the problem to a system of coupled differential equations for 
F{r) and G{r), 

dF fl + \ eBr\ ( Za\ ^ 

2^^-F+LE + me + \G = 0, (3.56) 



dr \ r 2/ \ r J 

dG f^~^2 ^Br\ ^ . ( -r^ Za\ 



dr ■ \ r +—)G+[E-m, + -)F = 0, (3.57) 

where a = e^ = 1/137 is the fine structure constant. Solving G{r) from the first equation 
in terms of F{r) and substituting into the second equation, we obtain the second-order 
ODE for F(r), 

„ /I 1 \ , r , , , , 2EZa + (l + ^)/ro 

F"+ i ]F'+ lE^-ml- eBil + 1 + ^ ^^^ 

\r t + Tq] y r 

fro + (/ + |)/ro (l + ir-iZaf (e5£^2 . ^ ^ q^ (g^g^ 



r + ro r^ 

- ATiTiKrincr i 

E+m 



where ro = ^f" . Applying the transformation, 



we obtain 



where 



F(r) = r« e-">i f(r), ? = ^/(i + ^T' - (Zaf, (3.59) 

/"+(Hi±i_^_,B,);. + £!r!±|!l±£i; = 0. (3.60) 

\ r r + ro / r\r + ro) 

C2 = E2-m2-e5(e + / + ^), 



ci = 2EZa + 



i^'-m^-eS(e + / + ^) 



^0, 



Co = 2EZaro + / + --6 (3.61) 

This ODE is of the form (II. 5p and has, from our general results in previous sections, 
polynomial solutions of degree n = 0, 1, 2, ■ ■ ■ , 

n 

j{r) = J](r - n), / = 1 forn = 0, (3.62) 

where rj are the roots of the above polynomial to be determined, provided that E, Z, B 
are given by 

E^ = ml + eB(n + l + i+A, 



n 



2EZa = e5 I ro + ^ 

1 ( "^ " \ 

2EZar^ = -n{n + 2^ - 1) + ^ - (/ + -) + e5 J^ rf + ro ^ r, , (3.63) 
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and the Bethe ansatz equations 

V ^ +^i±l 1 eBri = 0, i = l,2,---,n. (3.64) 

^ Tj - rj Ti Ti + ro 

Let us remark that our f l3.63p -f l3.64p differ from the corresponding equations (22)-(25) 
of ref. [25]. It would be interesting to establish the relation between the two sets of 
expressions. 

Note that /(r) = 1 is a solution of the ODE fl3.60p provided that E,Z,B obey the 
following relations: 

E^=ml + eB(l + ^+^j, (3.65) 

2EZa = eBro, 

2EZaro = ^-{l + -). (3.66) 

These relations can be obtained from fl3.63p by setting n = 0. Solving these relations we 
obtain 

^ mlJl + l + O 

ii + i + 0'' 

E = -!^ + -Jml + 2eB = -— ^^ -, (3.67) 

2 2^ ' 2(/ + l + 

where ^ is related to the parameter Z via the expression given in fl3.59p . It follows that 

2 (: / 1 _(_ pRr'^ 2 

F(r) = r« e-^^'r, G(r) = }- \ ^ r« e'^^'r (3.68) 

with ei? and E given by fl3.67p in terms of ^ (i.e. Z). As far as we know, fl3.68p gives the 
first exact solution to the planar Dirac electron system. However, this solution does not 
seem to be squarely integrable. 

3.5 Schrodinger equation of 0{N) invariant decatic anharmonic 
oscillator 

Consider the 0{N) invariant decatic anharmonic oscillator in A^ dimensions [26]. The 
Schrodinger equation is 

iv' + y(x) j ^(x) = E^(x) (3.69) 

with the potential l^(x) defined by 

N 

\/(x) = Aix^ + X, (x^)' + A3 (x^)' + A4 (x^)' + (x^)' , x^ = J2xl (3.70) 

1=1 
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Here without loss of generality we have normalized the potential such that the coefficient 
of (x^) equals to 1. In A^- dimensional spherical coordinates, the radial wave function 
R{r) satisfies 



ar'^ r ar r^ 



R{r) = 0. 



(3.7i: 



Applying the transformation R{r) = r 2 ■?/'(r) yields 

X/i- 1^ 



-r + 



+ 2{V{r)~E) 



4j = 0, 



(3.72) 



where fi = I + ^{N — 1). Applying the transformation 



%l) = r e 



/^.-«V-/3V-7V 



z = r , 



(3.73) 



where a,/3,7 > are parameters yet to be determined, yields the ODE for (l){z 
+ 



— 'fz — pz — a 



j 0' + ^ [(2«/3 - ^{N + 21 + 4)- 2A2) z^ 
+ {a^ - /3{N + 21 + 2)- 2Ai) z + 2E- a{N + 21)] = 0, (3.74) 



where a, 13,^ are given by 



A4 



" = ^' ^ = 7!' " V2 



hM) 



(3.75) 



The ODE (I3.74p is of the form (II. 6p . Applying our general results in previous section, we 
obtain that this equation has polynomial solutions of degree n = 0, 1, 2, ■ ■ ■ , 



= n(^ - Zi), = 1 for ra = 0, 



(3.76) 



i=l 



where Zi are the roots of the above polynomial to be determined, provided that 

2a/3 - -f{N + 21 + A)- 2A2 = 4:jn, 

n 

- f3{N + 21 + 2)- 2Ai = 47 ^ Zi + 4/3r 



a 



in. 



i=l 



E=^{An + N + 2l) + 2pJ2''i + ^^J2' 



(3.77) 
(3.78) 

(3.79) 



1=1 



and the Bethe ansatz equations. 



= l^i + I3zi + a 

7: 7 : 



l + N/2 



j¥=i 



Zi Zj 



i=l 



i = 1,2,- ■ ■ ,n. 



(3.80) 
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Note that </> = 1 is a solution to f l3.74p with energy 



E 



2^2 



A,- 



A^ 



{N + 21), 



(3.8i: 



where A3 and A4 obey the constraints, 



A4 Uz - x) ~ ^(^ + 2/ + 4) - 2A2 = 0, 

-^(A^ + 2/ + 2)-2Ai = 0, 



A3-- 
^ 4 



A/i \ A4 



V2' 



(3.82) 
(3.83) 



which can be obtained from (13. 771) - (13. 791) by setting n = 0. The real solutions of these 
constraint equations (so that the energy E is real) are given by 



A. 



3(A^ + 2/ + 2) V 2 



?)2 + (?)3 



1/3 



+ 



Q 



■2' ^3^ 



1/3 



A3 = ^ + 



\l V2{N + 2I + A + ^/2\2) 



(3.84) 



where 



P 



24A2 



9(A^ + 2/ + 2)2' 

32V2XI {N + 21 + 4 + y2A2)2 



27(A^ + 2/ + 2)3 N + 21 + 2 

We thus obtain the first exact ground state of the radial Schrodinger equation ( 13.7ip . 

1 



(3.85) 



R{r) = r' exp 



2^2 



A3-flr= 



A4 4 1 

r 



4^2 3^2 



(3.86) 



with A3 and A4 given by (I3.84p . 

For n = 1, we find that the real root to the Bethe ansatz equation (I3.80p is 



A4 f V 



|)=-(i)' 



1/3 



+ 1-2- 



|)^-(i)' 



1/3 



(3.87) 



where 



u 



A4 _ A3 

24 ~Y' 

432 2 ^ 2v^ 



{N + 21). 



(3.88) 
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The parameters A3 and A4 are determined from the equations (see fl3.77p and (13.781) ). 

A4(A3-^j -v^(Ar + 2/ + 8)-2A2 = 0, (3.89) 



The energy eigenvalue E is given by 

1 / A^\ 

E = ^ iXs - ^j {N + 21 + A) + V2X^z^ + 2V2zl 

and the radial wave function is 
R(r) = r' (r^ — Zij exp 



(3.90) 



(3.91) 



T.M) 



2 _ -^4 4 



(3.92) 



2V2 V 4 y 4^2 3^2 

with Zi given by (I3.87p . This gives the first excited state of the system. 

4 Concluding remarks 

The main results of this paper are theorem II. ![ corollary 11.31 and their applications to the 
exact solutions of the physical systems (examined in section [3]). 

As spelled out in theorem II. H we have found all polynomials Z{z) such that the ODE 
(II. ip has polynomial solutions S{z) of degree n with distinct roots zi,Z2,- ■ ■ ,Zn. We 
have also found a set of n algebraic equations which determine the roots Zi and thus the 
corresponding polynomial solutions S{z). If the coefficients of polynomials X{z) and Y{z) 
in (11.11) are algebraically dependent, i.e satisfy the relations in corollary 11.31 then ODE 
(II. ip allows an sl{2) algebraization. We have also found all the polynomials Z{z) and 
degree n polynomial solutions S{z) under these conditions. 

Having as many as 12 parameters, the ODE ( 11. ip is very general and contains, as 
special cases, most known second order differential equations which occur in physical, 
chemical and engineering applications in the literature (e.g. the relatively simple Heun 
equation and its various confluent equations). Thus our general results (theorem II. II and 
its corollaries) provide an unifled derivation of exact, closed form solutions for all such 
systems. As applications to theorem ll.il we have examined flve physical systems in section 
[31 which are described by the ODEs corresponding to (ll.2p -( lL6l) . respectively. We have 
shown that these systems are quasi-exactly solvable, i.e. the corresponding differential 
equations have polynomial solutions, if their parameters satisfy certain constraints (special 
cases of ( ]1.8p -( fLl0l) ). The quasi-exact solvabihty has enabled us to use theorem 11.11 to 
obtain the closed form expressions for the eigenvalues and eigenfunctions of these systems. 

Our results (theorem 11.11 and corollary II. 3p can be extended to second order ODE of 
the form (fHI]) with degX{z) < I, degY{z) < / - 1 and degZ{z) < / - 2 for / > 5 as 
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well as to higher order ODEs. Research on this as well as on applications of theorem 11.11 
and corollary 11.31 in various areas of science is in progress, and results will be reported 
elsewhere. 
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5 Appendix 

For completeness and convenience of applications, in this appendix we write down the 
explicit formulas obtained from applying the theorem II. II to the special cases ( ll.2p - (]1.6p . 

Corollary 5.1 The coefficients ci,co of Z{z) such that the Heun equation ( fi.^) has poly- 



nomial solution ( [i.?| j are given by 



Cl 



—n 



n 



-i+E 



a. 



s=l 



a. 



y^^Zj + n{n - 1)^4 



i=l 



s=l 



Co = - 2(n-l) + ^, 

+n[ai{d2 + rfs) + 0:2(^1 + d'i) + a-i{di + ^2)], 
where the roots zi,Z2,- ■ ■ ,Zn are determined by the Bethe ansatz equations, 



(5.1) 



(5.2) 



Ez ^^^ as 






Zi Ug 



0, z = 1, 2, ■ ■ ■ ,?T,. 



(5.3) 



Write Cl = af3. Then (15. ip is nothing but the so-called Fuchsian relation, a + (3 + 1 



X;,=i«s, where 



a 



-n and /3 = Ylis=i Ois + n — 1. 



Corollary 5.2 The coefficients C2,ci,co of Z{z) such that the generalized Heun equation 
( fi.gj) has polynomial solution ( [i.7| ) are 

C2 = -n \ 2J ^is + n-l 

4 



,s=l 

4 



Cl 



Co 



-[Y,^^s + '2{n-l)\Y,Zi + 



n 



4 



n 



[n 



l)^e, + P 



- J] /i. + 2(n-l) 5^^,2 + 2^ 



Z'lZj 



s=l 



i<j 



(5.4) 
(5.5) 
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2(n-l)^e, + P 



s=l 



^Zi + n{n - 1) ^ e^ej + Qn, 



(5.6) 



i=l 



s<t 



where 



P = /ii(e2 + 63 + 64) + /i2(ei + 63 + 64) +/i3(ei + 62 + 64) +Ai4(ei + 62 + 63), 

Q = Ail(6263 + 6264 + 6364) +/i2(6i63 + 6164 + 6364) 

+yU3(6l62 + 6164 + 6264) + /i4(6l62 + 6163 + 6263), 

and the roots zi, Z2,- ■ ■ , Zn are determined by the Bethe ansatz equations, 



fJ'S 



E^ + E ._ 



0, i = 1, 2, ■ ■ ■ ,n. 



(5.7) 



Corollary 5.3 The coefficients C2,Ci,Co of Z{z) such that the differential equation l^l.^ 
has polynomials solution ( [i.7| ) are given by 



C2 = 


= —nu, 


Cl = 


n 

-- -v^Zi-n 

s=l 


Co = 


n 
s=l 



(n-l) + ^z/, -z/^/. 



3 3 

s=l s=l 

3 3 

s=l s=l 



(5.8) 
(5.9) 



^Zi + n{n-l)^fs 



i=l 



s=l 



2{n-l) + ^Us-y^fs 

s=l s=l 

+ Hflf2 + /1/3 + /2/3) - l^l{f2 + /3) - 1^2{fl + /3) " ^^3(/l + /2)] ^i, (5.10) 

where the roots zi, 2:2, ■ ■ ■ , -2„ satisfy the Bethe ansatz equations, 



v^ 



jj^i ■> s=l 



+ 1/ = 0, i = 1, 2, ■ ■ ■ ,n. 



(5.11) 



Corollary 5.4 T/ie coefficients C2,Ci,Co of Z{z) such that the differential equation ( li.5|) 
/ias polynomials solution ([i.7| ) are 



C2 = -ncr, 

n 

Cl = -(T^2;i -n[fi;-or(5fi +5(2)] 



(5.12) 
(5.13) 



Co = -a^^f - [K-a(^i + ^2)]5]^2i 

-n(n - 1) - n [cTi + 0-2 + o-5'i5'2 - f^igi + 5'2)] , 
where the roots zi,Z2,- ■ ■ ,Zn obey the Bethe ansatz equations, 

n „ 

^ \ = h CT-Zj + K = 0, i = 1,2, ■ ■ ■,n. 



il^i 



Zi — Zj Zi — Qi Zi — g2 



(5.14) 



(5.15) 
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Corollary 5.5 The coefficients C2,Ci,Cq of Z{z) such that the differential equation U.6]) 
has polynomials solution ( [j.7| ) are given by 

C2 = —nX, (5.16) 

n 

ci = — A V^-2j — n(7 — A/i), (5-17) 

n n 

Co = -A5];z,2-(7-A/i)J]^,-n(5-7/i), (5.18) 

s=l i=l 

where the roots z\^ Z2., ■ ■ ■ ■, Zn are determined by the Bethe ansatz equations, 

n „ 

^^^ + ^-r + A;z2^7Zi + ,5 = 0, i = l,2,---,n. (5.19) 

. / . '^2 '^7 Zi n 
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